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SPECIFICATION 

BLOOD FLOW DYNAMIC ANALYZING APPARATUS 
AND METHOD, AND IMAGE DIAGNOSING APPARATUS 

5 Technical Field 

The present invention relates to a blood flow dynamic analyzing apparatus and 
method and an image diagnosing apparatus, and in particular, to a technique that is useful 
for using a blood flow dynamic analyzing apparatus or the like to execute image analysis 
on a tomogram provided by a computer tomogram diagnosing apparatus that can use an 
10 X ray or an electromagnetic wave to provide biological function information dependent 
on a temporal variation in a subject. 

Background Art 

A conventional blood flow dynamic analyzing apparatus uses, for example, an 
X-ray CT apparatus to carry out dynamic imaging. During the imaging, an 

15 iodine-based contrast medium is injected into a patient. The blood flow dynamic 
analyzing apparatus thus provides information on a temporal variation in the 
concentration of the injected contrast medium. A time-concentration curve can be 
obtained from the temporal variation. Then, by analyzing the curve of the 
time-concentration curve, it is possible to obtain biological function information such as 

20 blood flow dynamics. Typical algorithms for analyzing the blood flow dynamics 

include a gamma fitting method, a maximum slope method, and a deconvolution method. 

First, the gamma fitting method uses a gamma function to approximate the 
time-concentration curve. Then, blood flow information is calculated from a peak value 
of the approximate curve and an area under the curve. Next, the maximum slope 

25 method calculates a blood flow by dividing the maximum value of slope of the 

time-concentration curve for each tissue by the maximum value of a rise in a CT value in 
an arterial input function. 

However, the gamma fitting method and the maximum slope method require 
that about 8 to 10 ml/sec of contrast medium be injected into the subject. 

30 Disadvantageously, this is a heavy burden on the physical strength of the subject into 
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which the contrast medium is injected. Moreover, disadvantageously, the gamma 
fitting method enables qualitative evaluation but not quantitative evaluation. 

Thus, efforts were made to develop a blood flow analyzing method which can 
reduce the contrast rate and which enables quantitative evaluation. As a result, the 
5 decon volution method has been proposed (see the well-known document shown below). 
The deconvolution method subjects an arterial input function and a tissue residue 
function to deconvolution to generate an impulse residue function. Then, blood flow 
information is calculated from a peak value or an area under the curve of the impulse 
residue function generated. The deconvolution method advantageously enables 

10 examinations at a low contrast rate of about 3 to 5 ml/sec. Accordingly, the 

deconvolution method requires a contrast medium injection speed that is only about half 
that of the ganrnia fitting method or the maximum slope method. 

(Well-known document: L. Ostergarrd etc.: High Resolution Measurement of 
Cerebral Blood Flow using Intravascular Tracer Bolus Passages: 1996; Magnetic 

15 Resonance in Medicine Vol. 36: P. 715-725) 

However, the deconvolution method requires a plurality of integrations in the 
data conversion calculation that determines the impulse residue function from the arterial 
input function and tissue residue function. Consequently, the deconvolution method 
requires a time for the calculation. 

20 It is an object of the present invention to enable blood flow dynamic analysis to 

be carried out in a short time. 

DISCLOSURE OF THE INVENTION 
To accomplish the above object, the present invention is characterized by 
inputting a tomogram acquired by a tomograph, determining a time-concentration curve 
25 indicative of information on a temporal variation in each pixel of the input tomogram, 
determining a time-concentration curve for an inflow artery from the determined 
time-concentration curve for each pixel of the tomogram, generating an inverse filtering 
function from the determined time-concentration curve for the inflow artery, generating a 
transfer function for each pixel of the tomogram on the basis of the reverse function 
30 generated and the time-concentration curve for each pixel of the tomogram, using the 



BEST AVAILABLE COPY 



generated transfer function for each pixel of the tomogram to determine a blood flow 
dynamic analysis image. 

This makes it possible to, for example, generate a transfer function for each 
tissue by multiplying, on a frequency space, the inverse filter of the time-concentration 
5 curve for the inflow artery and a time-concentration curve for the tissue. It is thus 

possible to reduce the time for calculations compared to the conventional deconvolution 
method. 

Further, according to a desirable embodiment, the present invention is 
characterized by extracting a maximum connected pixel area from the tomogram, 
10 creating a mask image showing the extracted maximum connected pixel area, and using 
the mask image created to remove unwanted areas from the input tomogram which are 
different from the area of the mask image and which include room air, a bed, and bones. 

Thus, the areas unnecessary for the blood flow analysis are removed to reduce 
the amount of data to be analyzed. This makes it possible to reduce the time required 
15 by the above biological function information calculating means or the like to execute 
calculations. 

Furthermore, according to another desirable embodiment, the present invention 
is characterized by determining a time-concentration curve for an outflow vein from the 
determined time-concentration curve for each pixel of the tomogram and finding peak 

20 values of the time-concentration curve for the inflow artery and the time-concentration 
curve for the outflow vein, and by correction of a partial volume averaging (simply 
referred to as "PVA" below), effect in the time-concentration curve for the inflow artery 
such that the peak value of the time-concentration curve for the inflow artery matches the 
peak value of the time-concentration curve for the outflow vein. 

25 Specifically, if an X-ray CT apparatus is used to examine blood flow dynamics, 

the PVA effect reduces the CT value of a voxel containing the artery compared to its 
inherent value. This degrades the quantitativeness of the results of the analysis. 
Therefore, correction of the PVA effect enables a clear blood flow dynamic analysis. 

BRIEF DESCRIPTION OF THE DRAWINGS 
30 Figure 1 is a block diagram showing an example of the hardware configuration 

of a blood flow dynamic analyzing apparatus according to the present invention; 
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Figure 2 is a flowchart showing a procedure from inputting of a tomogram to the 
blood flow dynamic analyzing apparatus according to the present invention to displaying 
of a function image carried out by the blood flow dynamic analyzing apparatus; 

Figure 3 is a flowchart showing a method for removing unwanted areas; 
5 Figures 4(a) and 4(b) are diagrams showing a method for correcting the PVA 

effect according to the present embodiment; 

Figures 5(a) to 5(c) are diagrams showing another method for correcting the 
PVA effect according to the present embodiment; 

Figure 6 is a flowchart showing a procedure for calculating a vessel diameter 
10 and a peak value according to the present embodiment; 

Figure 7 is a conceptual drawing illustrating an inverse filter method according 
to the present embodiment; 

Figure 8 is a conceptual drawing illustrating the inverse filter method according 
to the present embodiment; 
15 Figure 9 is a flowchart showing a process procedure of an inverse filter 

calculation according to the present embodiment; 

Figure 10 is a diagram illustrating a method for calculating biological function 
information; 

Figures 11(a) to 11(c) are diagrams illustrating a method 1 for calculating a 
20 high-frequency control filter; 

Figure 12 is a diagram illustrating a method 2 for calculating a high-frequency 
control, filter; and . . - . - 

Figure 13 is a diagram illustrating a method 3 for calculating a high-frequency 
control filter. 

25 BEST MODE FOR CARRYING OUT THE INVENTION 

With reference to the accompanying drawings, a detailed description will be 
given of a preferred embodiment of a blood flow dynamic analyzing apparatus according 
to the present invention. 

Figure 1 is a block diagram showing an example of the hardware configuration 
30 of a blood flow dynamic analyzing apparatus according to the present invention. As 
shown in the Figure 1, the blood flow dynamic analyzing apparatus is composed of a 
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tomogram input device 1 to which a tomogram acquired mostly by a computer 
tomography imaging apparatus is input and a calculating device 2 that executes various 
calculations such as image analysis. 

The tomogram input device 1 and the calculating device 2 need not necessarily 
5 be separate from each other but may be integrated. 

The calculating device 2 is composed of an interface (I/F) 3 that loads a 
tomogram from the tomogram input device 1, a memory 4 in which tomograms and 
results of calculations are temporarily stored, a central processing unit (CPU) 5 that 
executes various calculations, a recording medium 6 such as a hard disk on which results 
10 of calculations and function images are recorded, an external input device 7 such as a 
mouse or a keyboard, and a display device 8 that displays results of processing. These 
components are connected together by a common bus 9. 

Figure 2 is a flowchart showing a procedure from inputting of a tomogram to the 
blood flow dynamic analyzing apparatus according to the present invention to displaying 
15 of a function image carried out by the blood flow dynamic analyzing apparatus. 

First, the tomogram input device 1 selects a tomogram of a subject on which a 
blood flow dynamic analysis is to be executed. The tomogram input device 1 then 
inputs the selected tomogram to the calculating device 2 (step 201). The input 
tomogram is temporarily stored in the memory 4 or saved to the recording medium 6. If 
20 the tomogram input device 1 and the calculating device 2 are integrated, the external 

input device 7 selects a tomogram of a subject on which a blood flow dynamic analysis is 

. - to be executed. The selected tomogram is read from the recording medium 6.and - 

temporarily stored in the memory 4. 

Subsequently, the CPU 5 removes areas in the tomogram which show room air, 
25 a bed, bones, and the like and which are unnecessary for analyzing biological function 
information (step 202). A method for removing unwanted areas will be described later. 

Then, the CPU 5 obtains a time-concentration curve indicative of information on 
a temporal variation, from the tomogram for each pixel (step 203). Subsequently, the 
CPU 5 obtains a first circulatory component from the time-concentration curve (step 
30 204). To obtain the first circulatory component, it is possible to use an arbitrary 
well-known algorithm such as gamma function fitting or extrapolation based on an 
exponential function. 

5 
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Subsequently, the CPU 5 determines an artery used for an artery input function 
(referred to as an inflow artery below) and a vein used for quantization such as correction 
of the PVA effect (referred to as an outflow vein) (step 205). An operator may 
manually specify an inflow artery and an outflow vein while viewing a tomogram. It is 
5 also possible to automatically determine an inflow artery and an outflow vein for an 
organ to be analyzed on the basis of a peak value and peak time of a time-concentration 
curve. 

Before the inflow artery and the outflow vein are automatically determined, 
central pixels of the inflow artery and outflow vein are first calculated. In order to 

10 automatically select the central pixels of the inflow artery and outflow vein, it is possible 
to utilize a difference between the maximum and minimum values of a 
time-concentration curve for each tissue (referred to as ATDC below) or a difference in 
peak time or peak value characteristic. For example, by selecting one pixel having a 
peak value equal to or larger than a predetermined threshold and having the earliest peak 

IS time from pixels having a ATDC equal to or larger than a certain threshold, it is possible 
to automatically select the central pixel of the inflow artery. Further, the central pixel 
of the outflow vein can be automatically selected by selecting one of pixels having a 
ATDC equal to or larger than a certain threshold, the one having a peak value equal to or 
larger than a predetermined threshold and having the latest peak time. Then, for pixels 

20 around the central pixels of the inflow artery and outflow vein, pixels having a peak 

value equal to or larger than the threshold are separated from those having a peak value 
smaller than the threshold.. Then, connected pixels including the central pixels are 
extracted from the pixels having a peak value equal to or larger than the threshold. As a 
result, an inflow artery area and an outflow vein area are determined. The central pixels 

25 alone may be determined to be an inflow artery and an outflow vein. 

Then, the CPU 5 corrects the PVA effect (step 206). A method for correcting 
the PVA effect will be described later. Subsequently, the CPU 5 executes an inverse 
filter calculation on each pixel to calculate biological function information (step 207). 
The reverse filter calculation will be described later. 

30 A function image is obtained by mapping biological function information 

analyzed for each pixel (step 208). The display device then displays the function image 
(step 209). Step 204 may be executed between steps 206 and 207. 



6 
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Now, a method for removing unwanted areas will be described. 

Figure 3 is a flowchart showing a method for removing unwanted areas. As 
shown in the Figure 3, the CPU 5 first binarizes the tomogram input in step 201 in Figure 
2 (step 301). The binarizing process is executed by, for example, replacing, in a 
5 tomogram, pixels having a value equal to or larger than the threshold with 1, while 

replacing pixels having a value smaller than the threshold with 0. The threshold may be 
such an arbitrary value as enables a living tissue to be analyzed to be separated from 
room air. For example, if the input tomogram is a CT image, the threshold may be set 
at about -200. 

10 Subsequently, the CPU 5 labels the image binarized in step 301 (step 302). 

The labeling process is executed by applying the same number (label) to all the 
connected pixels (connected components), while applying a different number to different 
connected components. If adjacent pixels both have a value of 1, they are connected 
together. Values applied as labels are, for example, consecutive numbers starting at 50. 

15 Then, the CPU 5 searches for the maximum connected components (step 303). 

Searching the maximum connected components is executed by scanning the entire 
labeled image, counting the number of pixels for each label value, and selecting 
connected components having a label value indicating the largest number of pixels. 
Subsequently, the CPU 5 extracts the maximum connected components (step 304). The 

20 maximum connected components are extracted by replacing, with 0, all the pixels except 
those selected in step 303 and having a label value indicative of the largest number of 
. . pixels. Thus, in the binarizing stage, pixel values for unwanted areas such as the bed - > 
and an intravenous drip infusion tube which may have their pixel values replaced with 1 
are replaced with 0. 

25 Then, the CPU 5 tracks the contour of the maximum connected components to 

extract the contour line of the outermost periphery of the maximum connected 
components (step 305). The contour line is tracked by horizontally scanning the image 
in which only the maximum connected components are left in step 304, starting with a 
pixel in an upper left comer of the image, using a pixel encountered first and having a 

30 nonzero label value as a start point to track the contour counterclockwise, and ending the 
tracking when the scan returns to the start point. A value different from the label value, 
for example, 1, is substituted for the pixels on the contour line. 
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Subsequently, the CPU 5 fills 1 into the inside of the contour line determined in 
step 305 to create a mask image (step 306). To fill the inside of the contour line with 1, 
it is possible to apply a seed fill algorithm that is a conventionally known process of 
filling the interior of a closed area with data using one point inside the closed area as a 
5 start point. In the present embodiment, the closed area is the contour line determined in 
step 305. Further, points inside the closed area are pixels with the label value for the 
maximum connected components determined in step 304. The inside of the contour line 
can be filled with data by detecting one pixel with this label value and using it as a start 
point to execute a seed fill process. 

10 Then, the CPU 5 multiplies the input image (tomogram) by the mask image 

created in step S306 to remove pixels located outside the mask area (step 307). 
Specifically, each pixel of the mask image is examined. When a zero pixel is detected, 
the minimum luminance is substituted for a pixel in the original image corresponding to 
the coordinates of the zero pixel. This operation is performed on all the pixels to 

15 remove the pixels located outside the mask area from the input image. 

Now, a method for correcting the PVA effect will be described. 
Figure 4 is a diagram showing a method for correcting the PVA effect according 
to the present embodiment. For example, in case that biological function information 
on the brain is analyzed, in general, the anterior or middle cerebral artery is used as an 

20 inflow artery, and the superior sagittal venous sinus is used as an outflow vein. 

According to the study of Lapin et al. (Journal of Computer Assisted Tomography 1993; 
Vol. 17: P. 108-144), the PVA effect reduces the signal value of a- voxel containing an 
artery, so that a time-concentration curve for the anterior or middle cerebral artery, which 
are an inflow arteries, has a value smaller than an inherent one. Further, according to 

25 the study of Lee et al. (American Journal of Neuroradiology 1999; Vol. 20: P. 63-73), a 
vessel of diameter at least 1.73 mm is not affected by the PVA effect, so that a 
time-concentration curve for the superior sagittal venous sinus, which is outflow vein, is 
not affected by the PVA effect. 



30 vein are generally as shown in Figure 4(a). A peak value Pa of the time-concentration 
curve for the inflow artery is smaller than that Pv for the time-concentration curve for the 
outflow vein by a value corresponding to a decrease in signal value attributed to the PVA 



Accordingly, the time-concentration curves for the inflow artery and outflow 




effect. Ideally, the concentration of a contrast medium in the vessel is constant 
regardless of its diameter. Consequently, in case of no influence of the PVA effect, the 
peak signal value is constant regardless of the vessel diameter. 

Therefore, when the peak value for the inflow artery is defined as Pa and the 
5 peak value for the outflow vein is defined as Pv, the PVA effect on the inflow artery can 
be corrected by multiplying the time-concentration curve for the inflow artery by Pv/Pa 
as shown in Figure 4(b). When a time-concentration curve for the inflow artery 
observed at a time t before the correction of the PVA effect is defined as AIFpva(t) and a 
time-concentration curve for the inflow artery observed at a time t after the correction of 
10 the PVA effect is defined as AIF(t), the correction may be made so that the formula 
shown below is established between these time-concentration curves. 



conditions as those for examinations to be carried out, the image can be used to correct 
the PVA effect. 

For example, in the clinical image shown in Figure 5(a), for a plurality of 
vessels 1 to 5 contained in the tomogram, vessel diameters Rl to R5 and peak values of 

20 time-concentration curves for the vessels are calculated. These values are shown in 

Figure 5(b). When these values are plotted in a graph, such a curve as shown in Figure 
5(c) is obtained. Description will be given later of a method for calculating the vessel 
diameter and the peak value. The thinner the vessel diameter is, the more pronounced 
the PVA effect is. The vessel is not subjected to the PVA effect when its diameter has 

25 at least a certain value. Accordingly, the curve in Figure 5(c) rises gradually until a 

diameter Rth is reached. Then, the peak value becomes flat for R equal to or larger than 
Rth. In this flat state, the peak value is defined as Pth for convenience. The vessel 
diameter of the inflow artery is defined as Ra. The peak value of the time-concentration 
curve for the inflow artery is defined as Pa. 

30 If Ra is smaller than Rth, the PVA effect on the inflow artery can be corrected 

by multiplying the time-concentration curve for the inflow artery by Pth/Pa. When a 



AIF(t) = AIFpva(t) — 



...(1) 



15 



Figure 5 is a diagram showing another method for correcting the PVA effect 
according to the present embodiment. 

If there is any past clinical image taken under the same imaging and contrast 
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time-concentration curve for the inflow artery observed at a time t before the correction 
of the PVA effect is defined as AIFpva(t) and a time-concentration curve for the inflow 
artery observed at a time t after the correction of the PVA effect is defined as AIF(t), the 
correction may be made so that the formula shown below is established between these 
5 time-concentration curves. 



Now, description will be given of a method for calculating the vessel diameter 
and the peak value. 

Figure 6 is a flowchart showing a procedure for calculating the vessel diameter 

10 and the peak value according to the present embodiment. First, the CPU 5 calculates 
the peak time and peak value of each time-concentration curve (step 601). Then, the 
CPU 5 calculates the central pixel of the vessel (step 602). The characteristics of the 
peak value may be utilized to calculate the central pixel. Selection of a pixel with the 
largest peak value enables the automatic selection of the central pixel of the vessel. 

15 Subsequently, for pixels located around the central pixel of the vessel, the CPU 5 
separates pixels having a peak value equal to or larger than a threshold from pixels 
having a peak value smaller than the threshold (step 603). Subsequently, the CPU 5 
extracts connected pixels including the central pixel from the pixels determined in step 
603 to have a peak value equal to or larger than the threshold. Thus, a vessel area is 

20 determined (step 604). 

- , ^ Then, the CPU 5 calculates the number- of pixels-in the vessel area (step 605). 
The CPU 5 then calculates the size of the vessel area from the calculated number of 
pixels (step 606). Since the area of one pixel has a fixed value in every apparatus, the 
size S of the vessel area is obtained by multiplying the area of one pixel by the number of 

25 pixels determined in step 605. 

Subsequently, the CPU 5 calculates the vessel diameter R (step 607). The 
shape of the vessel on the tomogram may approximated by a circle. Accordingly, the 
formula shown below holds true for the vessel diameter R and area S. The vessel 
diameter may be calculated in accordance with the formula. 



AIF(t) = AIFpva(t) 



... (2) 



30 




... (3) 
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Then, the CPU 5 calculates a time-concentration curve for the vessel (step 608). 
Since the time-concentration curve for each pixel was calculated in step 203, the 
time-concentration curve for the vessel is obtained by calculating the average of the 
time-concentration curves for all the pixels in the vessel area. Instead of taking the 
5 average, it is possible to carry out arbitrary weighting in accordance with the distance 
from the central pixel of the vessel calculated in step 602. 

Then, the CPU 5 calculates the peak value of the time-concentration curve for 
the vessel (step 609). It is possible to reverse the order of the processing from step 605 
to step 607 and the processing from step 608 to step 609. In the present embodiment, 
10 the vessel is extracted utilizing the characteristics of the peak value. However, the 

method for extracting the vessel is not limited to this. An arbitrary method may be used 
provided that it can accurately extract the vessel. Alternatively, the vessel may be 
extracted by manually specifying the vessel area while viewing the tomogram. 

Now, the concept of an inverse filter method will be described. 
IS Figure 7 is a conceptual drawing illustrating an inverse filter method according 

to the present embodiment. For a target tissue for a transfer function f, when an input 
function is defined as Cin(t) and an output function is defined as Cout(t), Cout(t) is 
obtained by subjecting Cin(t) and f to a convolution integral. This is expressed as 
follows. 

20 Cout(t) = Jcin(t) • f (t - f)df ... (4) 

0 

If a contrast medium is injected into the body through the artery, Gin(t) is 
replaced with the time-concentration curve (artery input function) AlF(t) for the inflow 
artery and Cout(t) is replaced with the time-concentration curve Cr(t) for the tissue. 
This is expressed as follows. 

25 Cr(t) = jAIF(t) • f (t - f)df . • - (5) 

0 

When Formula (5) is converted into a frequency space, the following formula is 

given. 

F{Cr(t)} = F{AIF(t)) • F{f(t)) ... (6) 
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In this formula, F{AIF(t)} denotes a Fourier transformation of AIF(t), F{Cr(t)) 
denotes a Fourier transformation of Cr(t), and F{f(t)} denotes a Fourier transformation of 
f(t). 

According to Formulae (5) and (6), the transfer function f can be generated by 
5 executing a reverse convolution on the time-concentration curve for the inflow artery and 
the time-concentration curve for the tissue. However, equivalent calculation results are 
obtained by multiplying the time-concentration curve for the tissue by the inverse filter 
of the time-concentration curve for the inflow artery instead of executing a reverse 
convolution. 

10 When the inverse filter of AIF(t) is defined as AIF \t), f(t) is expressed as 

shown below, on the basis of Formula (6). 



15 to the present embodiment. 

For blood flow dynamic analysis, in most cases, the contrast medium is injected 
into the body through the vein in order to avoid invading the subject. In this case, the 
contrast medium circulates through the body and then flows into the tissue to be analyzed. 
The input function is defined as Cin(t). A transfer function for a circulatory path is 
20 defined as fi. A transfer function for the target tissue is defined as f2. The output 
fiinction for the target tissue is defined as Cout(t). Then, Cout(t) is obtained by 

subjecting Cin(t) and fi to a convolution integral and subjecting the result and fa to a - - 

convolution integral. Thus, the following formula is given. 



f(t) = F*[F{AIF'(t)} . F{Cr(t)}] ... (7) 

In this formula, denotes a reverse Fourier transformation. 
Figure 8 is a conceptual drawing illustrating the inverse filter method according 




... (8) 



0 



25 



The convolution integral of Cin(t) and f i is replaced with the time-concentration 
curve (artery input function) AIF(t) for the inflow artery. Cout(t) is replaced with the 
time-concentration curve Cr(t) for the tissue. Thus, the following formula is given. 




0 



Formula (9) is equivalent to Formula (5), previously described. Accordingly, 
30 the transfer function fa is expressed as follows. 
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f2(t) = F'[F{AIF^(t)| . F{Cr(t)}] ... (10) 

Figure 9 is a flowchart showing a process procedure of an inverse filter 
calculation according to the present embodiment. As shown in Figure 9, the CPU 5 
executes a Fourier transformation on the time-concentration curve for the inflow artery 
5 (step 901). Then, the CPU 5 calculates an inverse filter from the Fourier-transformed 
time-concentration curve for the inflow artery (step 902). When the real part of the 
Fourier-transformed time-concentration curve for the inflow artery is defined as F{ AIFr} 
and its imaginary part is defined as F{ AIFr} , the inverse filter is expressed as follows. 

F{AIF-'(t)} = ^- 

[F{AIF,} + i{AIF,}] 

F{AIF^ } iF{AIF,} ^ 

[F{AIFJ^ + {AIF, }^] [F{Anv }^ + {AIF^ }^] 

10 Then, high-frequency components of the inverse filter are controlled (step 903). 

The high-frequency components of the inverse filter are controlled by filtering the 
inverse filter using a high-frequency control filter such as a Wiener filter or a 
Butterworth filter. Subsequently, the CPU 5 executes a Fourier transformation on the 
time-concentration curve for each tissue (step 904). The CPU 5 uses the inverse filter 

15 determined in step 902 to filter the Fourier-transformed time-concentration curve for the 
tissue (step 905). The CPU 5 thus generates a transfer function that executes an inverse 
Fourier transformation on the filtered Fourier-transformed time-concentration curve for 
the tissue (step 906). 

Then, the CPU 5 calculates biological function infomiation such as the blood 

20 flow dynamics from the transfer function determined in step 906 (step 907). A method 
for calculating biological function information will be described below. Steps 904 to 
907 are repeated a number of times equal to the number of pixels on the tomogram of the 
tissue to be analyzed. Additionally, it is possible to reverse the processing firom step 
901 to step 903 and the processing in step 904. 

25 Now, description will be given of a method for calculating biological function 

information from the transfer function. 

Figure 10 is a diagram illustrating a method for calculating biological function 
information. When the transfer function is defined as f(t), f(t) has, for example, such a 
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shape as shown in Figure 10. The blood flow (referred to as "BF" below) is detennined 
from the maximum value of f(t) as shown below. 
BF = f„ax ...(12) 

Further, the blood volume (referred to as BV below) is determined from the area 
under the curve of f(t) as shown below. 
BV = If(t)dt ...(13) 
The mean transit time (referred to as MIT below) is determined from the width 
of f(t) as shown below. 

rf(t)dt 



MTT 



= J£lI2£i ...(14) 



10 In this formula, f^ax denotes the maximum value of f(t). 

Now, description will be given of a method for determining a high-frequency 
control filter. In order to calculate a transfer function in accordance with Formula (9), it 
is necessary to execute an inverse convolution on the time-concentration curve for the 
artery and the time-concentration curve for the tissue. The inverse convolutions are 

15 roughly classified into methods such as singular value decomposition method which use 
a determinant and methods such as the inverse filtering method which uses a Fourier 
transformation. The method using a Fourier transformation can reduce the time 
required for calculations compared to the method using a determinant. With the inverse 
filtering method, the Fourier-transformed time-concentration curve for the tissue is 

20 filtered on a frequency space using the inverse filter calculated from the 

Fourier-transformed time-concentration curve for the artery. This calculation is a 
process for emphasizing the high frequency components. In general, biological signal 
components appear in a low frequency region, while noise components appear in a high 
frequency region. Accordingly, simple inverse filtering calculations unnecessarily 

25 emphasizes noise, making it difficult to accurately obtain biological function information. 
Thus, the high frequency components must be controlled for the inverse filter. 

It is assumed that the artery input function for the target tissue, that is, the 
time-concentration curve for the inflow artery, has been actually measured as shown in 
Figure 11a. The mean transit time (MTTAitery) of this time-concentration curve is 

30 determined as follows. 
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n 



2t-Ca(t) 



MTT. 



Artery 



n 



... (15) 



£ca(t) 



t=0 



In this formula, t denotes a time, n denotes the number of measuring points, and 



Ca(t) denotes a time-concentration curve. Such an ideally rectangular input function 
ideaI_Ca(t) as shown in Figure 1 1(b) is determined from the peak time (PT in Figure 1 1) 
5 of the time-concentration curve for the inflow artery and the mean transit time 

determined using Formula (16); the width of the ideally rectangular input function 
around the peak time is equal to the mean transit time. When Fourier-transformed 
ideal_Ca(t) is defined as F { ideaLCa(t) } , F { ideal_Ca(t) } is as shown in Figure 11 (c) . In 
this case, useful ecological information appears as low frequency components, while 

10 noise components that are insignificant for the analysis appear as high frequency 

components. Accordingly, for F { ideal_Ca(t) } , for example, as shown in Figure 1 1 (c), 
only the components from the angular frequency 0 to (27c)/k are left as they are, 
components from (2n)/k to (47i)/k are attenuated, and frequency components exceeding 
(47i)/k are cut for determining a shape of the filter. This high frequency control filter is 

15 given by the following formula. 



In this formula, co denotes an angular frequency, cos denotes a control start 
frequency, and coe denotes a control end frequency. When the width of ideal_Ca(t), that 
is, the mean transit time of Ca(t) is defined as k, cos = (2n)/k and ©e = (47i)/k. 

20 Accordingly, the shape of the high-frequency region control filter is determined from the 
time-concentration curve for the inflow artery. In the description of the present 
embodiment, for convenience, cos = (2n)/k and coe = (4n)/k, However, cos and coe are 
not limited to this. Arbitrary settings are possible; cos ^ (n7r)/k and coe = (2n)/k (m < n). 
In the present embodiment, cos and coe are calculated from the time-concentration curve 

25 for the inflow artery. However, the method for calculating cos and coe is not limited to 
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this. An arbitrary setting is possible using the maximum frequency as shown in Figure 
12. 

The high frequency control filter may be as shown below. 
filter((D) = filter2(a)) + k{filter2((D) - filterl(a))) 



filterl(a)) = 



filter2(o) 
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5 The filter expressed by Formula (17) suppresses the high frequency components 

while emphasizing the low frequency components. In this formula, k denotes a 
parameter indicating how the low frequency components are emphasized. If the 
following settings are made: cols = 0, ©le = 7i/k, ©2s = 0, (o2e = (37i)/k, and k = 0.75, 
such a high-frequency control filter as shown in Figure 13 is obtained. In the present 

10 embodiment, cols, cole, a)2s, and Q>2e are calculated from the time-concentration curve 
for the inflow artery. However, the calculating method is not limited to this. These 
values may be arbitrarily set on the basis of the maximum frequency. 

Further, the high-frequency control filter may display a graph showing the 
progress of measurements on an image display section provided in a computer tomogram 

15 diagnosing apparatus. In connection with this graph, it is possible to allow at least one 
of various parameters provided for the high-frequency control filter to be varied by using 
a well-known graphical user interface that uses a pointing device such as a mouse or a 
track ball to, for example, vary the shape of the displayed graph. The various 
parameters include the control frequency of the high-frequency control filter and the 

20 degree of the emphasis of the low frequency area. 

Thus, the original signal can be appropriately emphasized by varying the various 
parameters, for example, increasing the control frequency band of the high-frequency 
control filter or the degree of the emphasis of the low-frequency area. As a result, the 
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number of signal components increases to suppress the noise components. An optimum 
high-frequency control filter can thus be provided. This makes it possible to provide an 
image subjected to high-definition biological function information analysis for blood 
flow dynamics. 

As described above, according to the above embodiment, if biological function 
information analysis is to be executed on a tomogram dynamically acquired using a 
computer tomogram diagnosing apparatus, the analysis can be accomplished with a 
contrast rate lower than that of the gamma fitting method or maximum slope method and 
with a calculation time shorter than that of the decon volution method. The embodiment 
is also effective in preventing the degradation of quantitativeness attributed to the PVA 
effect. The embodiment is also effective in preventing an unnecessary increase in 
calculation time owing to unwanted areas and analysis precision without increasing a 
burden on the operator. 



17 



